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Abstract 

The whole class of minimally coupled and massive scalar fields which may be 
responsible for flattening of galactic rotation curves is found. An interesting rela- 
tion with a class of scalar-tensor theories able to isotropise anisotropic models of 
Universe is shown. The resulting metric is found and its stability and scalar field 
properties are tested with respect to the presence of a second scalar field or a small 
perturbation of the rotation velocity at galactic outer radii. 

Published in Astronomy and Astrophysics, 413, p 799, 2004 

1 Introduction 

One of the most fascinating cosmological problems is dark matter one: 99 percent 
of the Universe energy density would be hidden. A good indication of dark mat- 
ter presence is g iven by galactic rotation curves which disagree with Kepler laws 
lFreemaJl9 70 1 . Particularly, spiral galaxy rotation curves seem flattened at large 
radii. One possible explanation is that they are made of a luminous disk whose 
density exponentially decreased to adjust to a dark halo whose distribution evolves 
as y- 1 I Ru bin etalll997l) 

Dark matter nature is unknown today. From WMAP observations, we know that 
the matter of which we are made represents 4% of the Universe content, 23% is 
made of cold dark matter and 73% of dar k energy. The se exotic types of matter 
could be represented by scalar fields l Mat os et alJll999B . which are predicted by 
unification theories l Ellis e t alll99iih . Starting from this assumption, we will study 
how they could be responsible for the observed flattened rotation curves. Dark mat- 
ter is not the only wa y to explain them l Batt aner & Floridol20 00'). Hence, Milgrom 
iMilgr om 1983a b c), Sanders l Sanders 1990) and others have proposed to mod- 
ify new tonian theory (MOND) for galaxies outer radii and in Battane ~& Floridol 
1 1995), ad hoc magnetic field are considered. 

Hence, the physical framework of this paper will be the scalar tensor theories. 
Unification theories and particularly supersymmetry predict the existence of scalar 
fields, which thus deserve be taken into account in cosmology. Most of time , only 
their cosmo logical consequences are analysed: quintessence phenomenon iJohri 
1200 it Isahnill2000h . isotropisation (Fav 2001) or inflation for instance. However, 
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they could also be present at galactic scale. In this work, we are going to assume 
that the dynamics of galaxies at outer radii is described by a scalar tensor theory 
with a dust perfect fluid, neglecting the radiation. The scalar field <p will be mini- 
mally coupled, massive with a Brans-Dicke function representing its coupling with 
the metric. It is equivalent to consider that, at galactic scale, the gravitational func- 
tion is a constant but the potential U and the Brans-Dicke coupling function a> vary 
with <p. 

Concerning the geometrical framework, we will consider a spherical and static 
metric. These are reasonable assumptions since, generally, a galaxy has a rotation 
axis around which turn the stars with a velocity much smaller than light speed. We 
thus neglect dragging effects, justifying a static metric l Guzman & Matos 2000). 
We will be interested by galactic regions where rotation curves flatten and where 
most of the dark matter should be present, i.e. the galaxies outer radii. Indeed 
internal regions need few or not dark matter to explain their dynamics. 

Our goal will be to study the form of the metric and the scalar field proper- 
ties explaini ng why the galactic rotat ion curves are flattened. A similar work has 
been done in M atos et all feOOd) and lMatos & Guzmar] 1200 lh . In the first paper, 
a massive scalar tensor theory is studied with a fixed to but an unknown (J. After 
having found the metric compatible with flat rotation curves, it is shown that the 
only potential able to reproduce such a dynamics should have an exponential form, 
U = e**. In the second paper, with the same form for a>, two massive scalar fields 
and a perfect fluid are considered. One of the potentials is assumed to have an ex- 
ponential form and similar results are found. In the present paper, we are going to 
consider a single massive scalar field with a perfect fluid but both tu and U will be 
unknown functions of <p. Then we will look for the metric and relations between 
oj and U allowin g to get the observed flat rotation curves, thus generalising the 
results of lMatos et alJl200fT) to a larger class of scalar tensor theories. Moreover, 
we will test the stability of our results by considering a small perturbation of the 
rotation velocity or/and an additional scalar field. 

It is important to note that other types of rotation curve s exist, such as decre asing 
rotati on curves found by Casertano and van Gorkom ICasertano & van Go rkom 
Il99ll) . Moreover, it seems that bright compact galaxy rotation curves are slightly 
decreasing whereas low luminosity ones tend to be increasing, indicating that they 
have more dark matter. However in this work we will only take into account 
asymptotically flat rotation curves. Indeed, a large number of them seems to be 
well approximated by an Universal Rotation Curve IPersic etalli99ri whose for- 
mulation, adapted to spiral galaxies, tends to a constant at late times. It shows 
how rotation curves depend on galaxy luminosity: increasing or decreasing ro- 
tation curves would respectively correspond to low or high galactic luminosity 
but should tend to a constant at outer radii. It seems to be confirmed by Swaters 
lSwaterslfl999l who has examined a large number of dwarf galaxies and found 
that their rotation curves flattened over 2 disk scale lengths. MOND theories or 
the presence of electromagnetic fields can also predict this type of curves. Con- 
sequently, although all the rotation curves do not flatten, this type of behaviour is 
sufficiently observed or predicted to justify a particular attention. 

The plane of this paper is the following. In section 2, we look for metric form 
and scalar field properties allowing to get flattened rotation curves. In section 3, 
we discuss about these results. 
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2 Metric and scalar field mathematical proper- 
ties 



This section is divided in two parts. In the first one, we consider the presence of a 
single scalar field. We look for the metric and properties of the unknown functions 
at and U such that the rotation curves be flattened at galaxy outer radii. In the 
second one, we test our results stability by adding a second scalar field. We will 
use a static and spherical metric written as: 

ds 1 = -e^dt 2 + e 2A dr 2 + r 2 (dQ 2 + sin 2 8d(ft 2 ) (1) 

<p and A being some functions of r. 

2.1 With a single scalar field 

The action for a minimally coupled and massive scalar field with a perfect fluid is 
given by: 

S = j(R- ^i/r" -U+ ijz, m ) V^ 4 * (2) 

R is the Ricci scalar, if/ the scalar field, (j(if/) the Brans-Dicke coupling function 
and U(ifi) the potential. L„, is the Lagrangian describing a dust perfect fluid whose 
impulsion-energy tensor writes T" 13 = pifuP with u" the 4- velocity and p the den- 
sity of the dust fluid. We get the field equations and Klein-Gordon equation by 
varying the action with respect to the metric functions and scalar field: 



r- 2 [r(l-e- 2A )] 


= ^<A'V 2A + 1/2(7 +p 


(3) 


- r~ 2 (l - e~ 2A ) + 2r 


1,/ -2A W ,i2 -2A in,, 

2if/ 


(4) 


e- 2A ((f>" + (f>' 2 + tf>'A' - 


A'/2) = -^VV 2A -l/2t/ 


(5) 


e- 2K ip' 2 (ojij/-\ + 2e- 2A wif/-' 


[(2r-> - A' + (P')iy + f ] - = 


(6) 



A prime stands for a derivative with respect to r and a if/ indice, a derivative with 
respect to the scalar field. By subtracting equations I3l4t and by summing I4I5I . it 
comes: 

r(A' - <f)') - 1 + e 2A [l - \ j2r 2 (U +p)] = (7) 
r 2 A'(20' - 1) + 2r 2 ((f>" + <p' 2 ) + 4r<f>' +2- 2e 2A (l - r 2 U) = (8) 
Then, we derive (J and p as some functions of A and <f>: 

p = 1 l2e- 2A r 2 [-2 + 2e 2A + r(4 - r + 2rcf>')A' + 2r 2 ((f>' 2 + (/>")] (9) 

U = -l/2e~ 2 V 2 [2 - 2e 2A + r<p'(A + 2r<j>') + r 2 A'(2(f>' - 1) + 2r 2 (f>"] (10) 

We introduce these expressions in J5} to get at as a function of A, (ft and if/. Then 
putting the above forms of U, p and at in Klein-Gordon equation yields: 

A' [2(r - 2) + r(r - 12)^' - 2r 2 <ft' 2 ] - 20' [e 2A -3 + r 2 ((ft' 2 + (ft")] = (11) 

which is scalar field independent. Since we are interested by flat rotation curves, 
we assume that rotation velocity tends to a constant for large r. However, rota- 
tion curves as seen by an observer at infinity for a spherical symmetry, are given 



3 



hy Vrot = ■\/ r Sit,r/(2g„), as shown in iMatos et alJ ^2000^1 where a newtonian in- 
terpretation of this last expression is given. It implies that <p' — > V 2 ol r~ [ and then 
e 2 * — > r 2v ' 2 <". To simplify our results below, we define the following constants: 



C\ 


= 2{\ + Vl,-2V%,) 


c 2 


= C-l 


C3 


= -2(2 + 6V 2 , f + 


C 4 


= 


C5 


= -2l/ 2 , f (-3 - + Vtj(2 + 6V; ol + O 


Ci 


= _ 3 _ + v 4 


C7 


= _4 _ rot/ 2 - 2V 4 

^ L ^ v ml ^ v ml 


c» 


= _2 - 41/ 2 - 4V 4 



Introducing </>' in 1 1 H and integrating, we find for A: 



(12) 



Ao is a positive integration constant. This last expression is only physically mean- 
ing for large r where it tends to the constant [Aoc^ 5 - lj as l/r. This constant 
must be positive otherwise e A is not defined for large r and moreover, numerical in- 
tegrations seem to show that A diverges for a finite value of this coordinate. Then, 
from 0, ED and (EJ, we calculate that asymptotically: 



•1 + A, 



(ci + c 2 r)(c 3 r 1 + c 4 ) C5 



0- 



C3 + c 4 r 



U = 2{c 2 - 



(2c 4 - 3)A [c 8 + r(c 4 - 1)] (c 4 + ca^'f 



(13) 



(14) 



c 3 + c 4 r 

For large r, p and U vanish as r" 1 . This asymptotical behaviour for the perfect fluid 
energy density is the same as the nonsingular isothermal profile, one of the most 
frequent halos. The metric descri bing the galaxies ou ter radii where the rotation 
curves flatten is thus the same as in lMatos et alJ ^OOO) whatever a>. Considering a 
perturbation S(r) of V rot does not modify these results as long as rS' — » 0. 
From the form of the metric and since we have left a> undetermined, we can get for 
large r a relation between ll> and U such that the rotation curves be flattened. By 
summing J5J and J4j and taking into account asymptotical behaviours for p and U, 
we find the following three limits: 



, 12 , -1 
toif/ if/ 

U 

U' = U*ifr' 



-2Uir~ 



(15) 
(16) 
(17) 



c 4 - 3 + - — -p^T — and C/j are some constants. We use <16> and 1171 for 

respectively introduce U and replace if/' in <15> . Then, considering g and k two 
functions of r and rewriting 1 1 51 and 1171 as respectively a>i[/' 2 if/~ l — > g(r) and 
Uftj/ — * k(r), we have 

ajk 1 



4 -4 
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Using Jl 61 to replace r 4 and introduce U, it comes 



UU 2 TJ 2 8(r) -4 
^ l W) 2 



Since here = At 2 r 2 and £(r) = —2U\r 3 , we find: 

& * (18) 



For a given form of L'(i^), 1 16i defines a unique form for r(<f). U(<fs) and r(ip) 
may then be introduced in 1 17t . defining a unique form for if/ (<//). Then, r(ip) and 
(/''((/') may be introduced in <15> . defining a unique form for w(ifi). Consequently, 
for a given U(ip), < 1 51 161 defined a unique to(tf/). The same remark applies to i 1 8i . 
Consequently, for a given U, the system 1151171 or 1 1 8i uniquely define to and thus 
the class of scalar tensor-theories responsible for the rotation curves flattening for 
given potential or Brans-Dicke coupling function. However, only Jl 51 171 uniquely 
define if/(r). 

We have shown above that p asymptotically behaves as r~ 2 . Examining the equa- 
tion J3j, we note that the scalar field energy density p^ shall be written as ^ip' 2 e~ 2A + 
1/2U. Knowing the asymptotical limits of each of these terms, we deduce that for 
large r, p oc p^: the scalar field is quintessent. 

The limit ITSl is doubly important. Firstly, in iFavlteOOlLfcoOStlFav & Luminetl 



2003) it has been shown that a necessary condition for isotropisation of Bianchi 



ipUf: 

models was —r* — > i l , I being a constant in a close interval depending on the pres- 
ence of curvature and perfect fluid. Consequently, galactic scalar field properties 
for large r could match a cosmological scalar field present in the entire Universe 
which would allow for its isotropisation. Secondly, specifying one of the unknown 
functions to or U, J 1 81 allow determining in a unique way the other one: this limit 
gives a necessary and sufficient relation between these two functions such that the 
galactic rotation cu rves for outer radii be flattened. It thus generalises the work of 
iMatos etalJfeOOOh for which to was a known function of the scalar field leading to 
an exponential potential. 

In the following section, we examine the stability of these results with respect to 
the presence of a second scalar field. 

2.2 With 2 scalar fields 

When two massive scalar fields are present, the action may take the following 
form: 

S = f(R- f - ^ 2 .X -U+ -£-L m ) V=irf 4 x (19) 

The i/r, are two scalar fields such that coi = coi(if/ t ), co 2 = to 2 {<l> 2 ) ar >d U = 
U(!/fi,ij/2). This form of the action is not the most general one but allows testing 
the results of the previous section. The field equations are: 

r- 2 [r(l - e- 2 *)] = if ^V 2A + if ^V 2A + 1/21/ +p (20) 
Wi 2if/ 2 

_ r -2 (1 _ e -2K) + 2r-y e -2A = if ^V 2A + if tfe^ - 1/21/ (21) 

e -^ (f + + ¥A > _ A7 2) = -if ^f e -2A - if tf e -^ - 1/2U (22) 

2(4-i 2if/ 2 ' 

e-^imfl 1 )^ + 2e- 2K co^\ l [(2r" ] - A' + + tf] -11^=0 (23) 
e- 2 ^i{co 2 r 2 \ 2 + 2e- 2k co 2 r 2 l \{2r~' - A' + <p')i(, 2 + ift] -U f2 =0 (24) 



5 



Making the same calculus as in section lzTl we get an equation similar to 1 1 1 i . i.e. 
independent on the scalar fields: 

4 - 4e 2A + 4rV 3 + 2r 3 A' 2 (20' - 1) + r*A" - 2r<£'(4 - 2e 2A + r A")- 

4r 2 </>" + 2rM [6 - 2r - (r - 16)r<p' + 4r 2 0' 2 + r 2 (f>"] - 2?$" = (25) 

Equations forp and U are the same as j5J and i 10i . Equation I25i does not depend 
on U, a>[, a>2 or the scalar fields forms. Moreover, we always have — > V 2 ot r~ [ 
which asymptotically characterises a flat rotation curve. The solution for A is- 
sued from equation <25> is thus independent on the scalar fields and the unknown 
functions w, and U. It will be always the same, whatever U, a>\, a>2 and fa- In 
particular, if we consider the special case where one of the scalar fields is negligi- 
ble, one have to recover the same asymptotical form for A as when only one scalar 
field is present. Hence, the asymptotical solution for equation J25> should be the 
same as for Jl U : when 2 scalar fields are considere d, A tends to a constant a s r -1 
and A' vanishes as r~ 2 . This is in agreement with Matos & Guzman! 1200 ll) and 
implies that U and p also vanish as r~ 2 . These results are the same if we consider 
a perturbation 6 for the rotation velocity as long as 6'r and 6"r 2 are asymptotically 
vanishing. 

Anew, we find the following limits allowing to determine if a relation exists be- 
tween ll>2 and U when the rotation curves flatten: 

w^Vr' + W21A2V2 1 -» 2rV 2 (26) 
U -> U,r- 2 (27) 
V = Uftifi + UfrPi -» -2U ir - 3 (28) 

Let us put that u>jif/' 2 iftj l — > g, and t/^.^J — > fc;. Then, gi + g2 — > r~ 2 , ki + fc - > '" 3 
and we have 

UjU 2 2 gi _ 4 

It implies that only one of the g,(or fe) have to tend to r~ 2 (respectively r~ 3 ), the 
second one varying as or slower than this last limit. Moreover, equations J23I24I 
may be written as: 

( ^;V 2A y + 2^;V 2A (^ + a,') - u lh ^ = 

fa fa r 

Hence, since A — > const and <p' — > r , if g, < r~ 2 , ki have to vary slower than r~ 3 . 
We thus distinguish 2 possible behaviours for g, and £, : 

• case 1 : g, — > r~ 2 and A:, — » r~ 3 

As previously, it comes that \ and ^ tend to some constants. 

• case 2: gi — » r~ 2 , g 2 « r~ 2 , k\ — » r~ 3 et fc 2 << r 3 

Consequently, the dynamical effects of 1^2 are asymptotically negligible in 
the field equations and the metric functions dynamics does not depend on it. 
We find that m , tends to a non vanishing constant and - diverges or 
vanishes. 

We conclude that the scalar fields ij/, which are not asymptotically negligible are 

such that — jTr tend to some non vanishing constants. Hence, the results of the 
previous section are not modified by the introduction of a second scalar field. 



6 



Anew, we note that the scalar fields energy density which shall be ■^>f'f e 2A + 
^4 l 2 e ~ %K + 1/21/ asymptotically tends to r~ 2 and behaves as the one of the dust 
fluid. It means that the two (or the dominant) scalar fields are (respectively is) 
quintessent. 

3 Discussion 

In this work, we have looked for the characteristics of the metric functions and 
scalar field such that galactic rotation curves flatten for outer radii. For the metric 
we have got the following result: 

Let us consider a minimally coupled and massive scalar field defined by a po- 
tential U and a Brans-Dicke coupling function a) with a spherically static metric. 
When, at outer radii, the galactic rotation curves flatten, the metric is asymptoti- 
cally defined by ds 2 = -r 2v >-«'dt 2 + A.\dr 2 + r 2 {d0 2 + sin 2 9dif> 2 ), Aj being a constant, 
whatever the form of the potential or Brans-Dicke coupling function. 

This result is stable related to a small perturbation 6 of the rotation velocity such 
that rd' — > or if we consider a second scalar field if> as defined by the Lagrangian 
1191 . In this last case the perturbation must be such that S'r and 6"r 2 vanish for 
large r. It is in accorda nce with the asymptotical form of the metric found in 
iMatos & Guzmart EoOll) where an exponential potential was found to explain the 
flattening of the rotation curves. Here, this property is generalized to any functions 
a> and potential U with the following characteristics: 

Let us consider a minimally coupled and massive scalar field defined by a potential 
U and a Brans-Dicke coupling function oj with a spherically static metric. When, 
at outer radii, the galactic rotation curves flatten, the energy densities of the perfect 
fluid and scalar field vanish as r~ 2 : the scalar field is asymptotically quintessent. 

Moreover, a> and U are asymptotically related by the relation ^-S — > I 2 , I 2 being 
a constant, and U vanishes as r~ 2 . 

Hence the energy density of the scalar field shows that we are using an isothermal 
profile which decays in r~ 2 and fit the flat galactic rotation curves quite well and not 
a Navarro-Frenk- White l Navarro et al. Il997t) profile where the density goes like 
r~ 3 in the asymptotic region and whose corresponding metric and rotation velocity 
has been recently determined in Matos & Nunez 1 2003). When 2 scalar fields are 
present, again their energy density behaves as the one of the perfect fluid. Con- 
sequently, at least one of them is quintessent. The presence of q uintessent scalar 
fields in spiral galaxies have been examined in M atos & Guzmanl (2000) where the 
agreement with the observed rotation curves is shown. Moreover, the scalar fields 

which are not negligible are such that — tends to a constant, leaving the above 
last result unchanged for these scalar fields. 

Let us make some remarks on the quantity Firstly, any necessary con- 
dition expressed with the unknown functions 1 of a scalar tensor theory such that 
the metric converge toward a determined form must be invariant with respect to 
a scalar field transformation. Indeed, considering F{U, a>), a necessary condition 
such that ds 2 always tends to a determined form, since a transformation if/ = TC¥) 
of the scalar field keeps the metric, it must be the same for the necessary condition 

i.e. oi and U for the present case. 
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F(U,oS). One can easily check it is the case for Particularly there is a scalar 
field transformation which allows to rewrite the metric under the form 

/l 671" 
(R - lyiw - U + —L m ) V=^x 

corresponding to the one of M atos etalj fooOO) and leading to their results. How- 
ever, for most of u> functions, this transformation is not defined or analyticaly 
workable and thus the results of Mat os et alJ l2000l) can not be arbitrarily applied 
to any forms of couple (cj, U) whereas it could be important to keep both oj and 
U as depending on the scalar field. As instance, if the potential vanishes, thus 
mimicing a vanishing cosmological constant, compatibility of the theory with PPN 
parameters requires that oj — » 00 and oj^oj^ — > iNordvedtl 1 968tlWagoneJ 1970). 
which does not fit with a constant oj got after field redefinition. Indeed, it recovers 
the general problem of finding a metric whose dynamics is agreed with the ob- 
servations and whose potential is in accordance with, as instance, particle physics 
predictions for the form of the potential: for this, it is necessary to keep the free- 
dom of choosing a form for o>((f>). Keeping oj as an undetermined function of the 
scalar field thus allows finding the set of all scalar tensor theories able to produce 
flat rotation curves for any forms of oj and if, even when the above scalar field 
redefinition can not be anatically performed. 

i/rUl 

The part of the results concerning the convergence of to a constant seems 
strangely correlated to isotropisation of homogeneous models which also needs 
this condition l Favll200lL 12003: Fav & Luminet 2003). It shows that the proper- 
ties of a galactic scalar field allowing the flattening of rotation curves could match 
those of a cosmological scalar field favouring Universe isotropisation. 

Starting from the form of the potential, the property — > l~ allows recover- 
ing the Brans-Dicke coupling function such that the rotation curves could flatten 
and vice-versa. Let us examine some of the most studied potentials. Hence, if we 
consider an exponential potential U = , we find that to = k 2 t~ 2 \$i, in accordance 
with the results of Matos et al. 12000) as a particular case of the class of theories 
we have found, and if/ cc lnr. If we take U = ifr, the Brans-Dicke coupling func- 
tion should be oj = k 2 {~ 2 i[/~ l and 1// oc r~ 2k . Thus, we see that considering two 
unknow n functions oj and U instead of a single one leads to an important generali- 
sation of Matos et al. ( 2000). Indeed, assuming asymptotically flat ro tation curves 
fix the behaviour of one of the metric function, i.e. (f>. Consequently, in Mat os et all 
(2000), since oj is chosen, the potential is uniquely determined whereas in this pa- 
per we can find it depends on oj, thus explaining the asymptotical relation between 
oj and U . It follows that an exponential potential is not the only one allowing to get 

flat rotation curves but a whole class of scalar tensor theories such that — > £ 
leads to this property. We remark that a cosmological constant can not explain 
why the rotation curves flatten since the potential must evolve as r~ 2 . This is in 
agreement with the theories which try to solve the cosmological constant problem 
by considering it as a variable function rather than a true constant. 

The interpretation of scalar fields properties found in this work may be done 
at cosmological or galactic scales. Phenomena giving birth to scalar fields prob- 
ably have a cosmological nature and are related to particle physics theories. To 
our knowledge it does not exist any way to generate them by galactic process. 
Their properties based on rotation curves observations are coherent with their 
usual cosmological picture i.e. their quintessent nature and their role in the Uni- 
verse isotropisation. However, if we consider an asymptotically increasing rotation 
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curve, it could mean that tp' increased faster than r -1 . It would come from equa- 
tions <10t and {5J that the quantity U jp could diverge. If such rotation curves were 
observed (it is t he case but some doubts subsist on the fact that they could asymp- 
totically flatten IPersic et alJll99d) 1. it would mean that scalar fields, at least at 
galactic scale, would not be quintessent. Then the quintessence properties would 
only be valid for some types of galaxies. In Casertano & van Gorkom ( 1991), de- 
creasing, increasing or flattening rotation curves are studied. Observations seem 
to show that the first ones appear for bright galaxies and the second ones for faint 
galaxies. A possible interpretation would be that flat rotation curves would be the 
outcome of an equivalent mixing between luminous and dark matter, the decreas- 
ing or increasing rotation curves resulting of a respectively luminous matter or 
dark matter domination. This flat rotation curves interpretation is coherent with 
the scalar field quintessence property found in this paper. 

To conclude, this work generalises those of M atos et alJfeOOd) and lMatos & Guzmanl 
1200 ll) . The metric got in these papers and the quintessent nature of the scalar fields 
have been generalised for any form of a>. A relation between U and uj has been 
found and allows getting easily one of these quantities from the other. It selects 
the class of scalar tensor theories which could be in agreement with flat rotation 
curves and shows that an exponential potential is not the only one able to produce 
such a dynamics. Moreover, we have remarked that this class is also in agreement 
with Universe isotropisation. The stability of these results with respect to a small 
perturbation of rotation velocity or the presence of a second scalar field have been 
tested. A next step would be to consider a non minimally coupled scalar tensor 
theory, i.e. with a variable gravitational constant or/and a magnetic field that could 
play a fundamental role at a galactic scale. 
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